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SUMMARY 

Specialty methods are presented for the computational simulation of spe- 
cific composite behavior. These methods encompass all aspects of composite 
mechanics, Impact, progressive fracture, and component specific simulation. 

Some of these methods are structured to computationally simulate, In parallel, 
the composite behavior and history from the Initial fabrication through several 
missions and even to fracture. Select methods and typical results obtained 
from such simulations are described In detail In order to demonstrate the 
effectiveness of computationally simulating (1) complex composite structural 
behavior In general and (2) specific aerospace propulsion structural components 
In particular. 


INTRODUCTION 

Fiber composites have distinct advantages over conventional materials 
which make composites highly desirable candidates for application In aerospace 
propulsion structures. Some of these advantages Include the well known high 
strength/density and modulus/density ratios. Other advantages Include high 
fatigue resistance, resistance to defects, high damping, and "zero" thermal 
expansion. Practically all of these advantages derive from the Inherent 
anisotropic and heterogeneous structure of fiber composites. This Inherent 
structure which provides composites with their distinct advantages also sub- 
stantially Increases the complexity required to formally describe their struc- 
tural behavior. 

Over the past 15 yr NASA Lewis Research Center has been developing spe- 
cialty methods for the analysis/design of composite structures. By virtue of 
their complexity, these methods are computationally Intensive and, by neces- 
sity, are In the form of Independent computer codes or are embedded as modules 
In more general structural analysis codes (refs. 1 to 3). Collectively, the 
development of these various specialty methods constitute a new discipline 
Identified as computational composite mechanics. These methods have been 
structured to simulate the history of composite thermal, structural and mate- 
rial degradation behavior beginning with the fabrication processes through 
several flight missions. The objective of this paper Is to briefly describe 
these methods and to Illustrate their application to general composite behavior 
as well as to select aerospace propulsion structures. 

Briefly, these methods encompass all aspects of composite mechanics, 
Impact, progressive fracture, and component specific composite structural ana- 
lysis. The select applications Included In the paper were chosen to Illustrate 
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(1) the computational simulation of composite structural behavior of Importance 
to aerospace propulsion structures and (2) the structural behavior of specific 
composite structural components In aerospace propulsion service environments. 
The applications are presented with a brief description of: (1) the Importance 

of the behavior to be simulated, (2) the simulation procedure, and (3) typical 
results obtained and their respective significance. 

Interply Laver Progressive Degradation Effects 
on Composite Structural Response 

Fiber reinforced composites are optimally utilized when they are required 
to resist and/or transfer only In-plane loads. Fibers are oriented In dif- 
ferent layers (plies) to resist these In-plane loads either In tension or com- 
pression. The different plies are held together by the Interply (Interlaminar) 
layers (matrix) which provide the composite with the structural Integrity 
required to resist the In-plane loads. Any degradation of the Interply layer 
will, therefore, affect the composite structural Integrity. As a result, 
Interply layer degradation effects on fiber composite structural Integrity Is 
receiving considerable research attention. Interply degradation Is generally 
either In the form of delamlnatlons or progressive weakening of the Interply 
layer. In either case, the angleplled laminate will most likely behave like a 
stack of Individual layers (beams) Instead of a composite. Compared to the 
undelamlnated composite, Individual delaminated layer behavior dramatically 
degrades out-of-plane structural response such as flexural deflections (bend- 
ing), buckling, vibration, and Impact. For example: (1) the flexural deflec- 

tion Increases as the square of the number of delaminated layers; (2) the 
flexural stress Increases linearly with the number of delaminated layers; (3) 
the buckling load decreases as the square of the number of delaminated layers; 
and (4) the frequency decreases linearly with the number of delaminated layers. 

Computational simulation (ref. 4) was used to determine and assess the 
effects of Interply layer progressive weaklng (degradation) on various struc- 
tural responses of composite beams. The structural responses of Interest 
Include: (1) bending, (2) buckling, (3) free vibrations, (4) periodic excita- 
tion, and (5) Impact. The computational simulation consisted of composite 
mechanics In conjunction with finite element analysis. The Interply layer 
degradation effects on the various structural responses were determined and 
assessed as a function of the Interply layer modulus, varying from 1 m ps 1 down 
to 1 ks 1 , and even lower for some limiting cases. The results obtained show 
that Interply layer degradation generally has negligible effect on composite 
structural Integrity unless the Interply layer modulus degrades to about 10 ks 1 
or less. 

Computational simulation results obtained for the Interply (adhesive or 
matrix) layer degradation effects on the forced vibration response of a simply 
supported beam are presented In figure 1. The third mode vibration amplitude 
Is shown without damping and with damping as a function of Interply layer mod- 
ulus and shear stress. These results Illustrate that the Interply layer degra 
datlon effects on complex composite structural response can be quantified by 
suitable computational simulation. 

Dynamic Structural Response of Unnotched and Notched Fiber Composites 

Unnotched and notched flexural specimen testing continues to be popular 
In the composites community for characterizing and/or qualifying fiber compos 
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1 tes for structural applications. These experimental techniques are popular 
because of:. (1) simplicity, (2) adaptability to adverse environments, (3) 
availability of simple equations for data reduction, and (4) availability of 
written ASTM standard testing procedures. In addition, flexural specimen test- 
ing forces the material to respond like a structure by simultaneously subject- 
ing It to tensile, compressive, and shear stresses. Furthermore, notched 
flexural specimen (Charpy Impact test specimen) testing enables easy determin- 
ation of fracture toughness and Impact resistance. 

Fracture, in general (be It static, quasistatic, or dynamic) Is a dynamic 
event and, as such, Is a very complex process. The fracture progression In a 
flexural specimen Is controlled by the local dynamic stress field which Is also 
complex. The local dynamic stress field Is characterized by stress waves which 
are normal, shear, flexural, and surface waves. Each of these dynamic stresses 
can Initiate failure (or a defect) at a point and/or propagate this defect to 
fracture. For a complete understanding of these phenomena, a detailed descrip- 
tion of the local dynamic stress field Is a prerequisite. 

The detailed description/quantification of the local dynamic stress field 
can be determined using computational simulation. The computational simulation 
(ref. 5) consists of composite mechanics, and finite element transient ana- 
lysis. The procedure Is demonstrated by applying It to the Impact response 
(displacements, velocities, stresses, and wave propagations) In Charpy-type 
composite beams. 

For the geometry (fig. 2) and type of materials studied, the static ana- 
lysis results gave close conservative estimates for the dynamic stresses. 

Other Interesting results from the study are: (1) the Impact-Induced effects 

propagate In S-GLASS/EPOXY specimens sooner than KEVLAR/EPOXY or T-300/EP0XY 
specimens, (2) the vibration mode shapes are coupled, exhibiting flexural, 
stretching and shearing motions, (3) the specimen responds for the Impact load 
duration only, (4) severe stress magnitudes are developed at the notch tip 
region leading to fracture Initiated by combined transverse normal and shear 
stresses and followed by fiber breakage caused by large flexural tensile 
stresses In this region, lyplcal results are shown In figure 3 for the dis- 
placement propagation In beams made from three different composite materials. 
These results demonstrate the effectiveness of computational simulation to 
determine the dynamic response of composite structural components. 

Interlaminar Fracture Toughness In Fiber Composites 

Interlaminar delamlnatlon of composites Is a type of fracture mode which 
needs to be carefully examined and properly accounted for In the design of com- 
posite structures. Regions prone to delamlnatlons Include free edges, loca- 
tions of stress concentration, joints, damaged areas, and defects resulting 
from the fabrication process. 

One way to properly account for Interlaminar delamlnatlon In a design Is 
to determine Interlaminar fracture toughness parameters and then evaluate 
stress states which are likely to Induce Interlaminar fracture. Several test 
methods to determine fracture toughness have been proposed and are currently 
being used. However, these test methods may not discriminate the participating 
and respective magnitudes of the different fracture modes. An alternate 
approach Is a suitable computational simulation. 
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Recent exploratory research at NASA Lewis (ref. 6) resulted In the deve- 
lopment of a computational method (procedure) for simulating Interlaminar frac- 
ture toughness In fiber composites. This computational procedure Involves 
three-dimensional finite element analysis, Including modeling of Interply 
layers. In conjunction with the strain energy release rate concept and compos- 
ite micromechanics. The procedure Is suitable for determining global and local 
Interlaminar fracture toughness parameters as well as critical values of these 
parameters (fig. 4). Predicted results using this procedure are In good agree- 
ment with limited available data from end-notch flexure and/or mixed-mode frac- 
ture test specimens (fig. 5). 

During the development of the procedure for the computational simulation 
of fracture toughness It was observed that (1) fracture toughness Is a struc- 
tural parameter and not a material characteristic as Is normally considered, 
and (2) critical fracture toughness parameters for rapid defect growth In a 
structural component can be determined by computational simulation. Both of 
these demonstrate, In part, that computational simulation Is not only applica- 
ble to quantify complex structural behavior but It Is also applicable to Iden- 
tify new phenomena and/or relationships. 

Progressive Fracture In Fiber Composites 

Determination of progressive fracture In composite laminates Is funda- 
mental and generic to developing the methodology for quantifying composite 
structures durability and reliability. Research activities on progressive 
composite fracture at NASA Lewis during the last 10 yr have culminated In the 
development of the C0DS1RAN computer code (refs. 7 to 9). CODSTRAN (Composite 
Durability Structural Analysis) has been specifically developed for the compu- 
tational simulation of progressive fracture In fiber composites. 

CODSTRAN Is a modular program that does quantitative calculations to pre- 
dict defect growth and progressive fracture In composite structural components. 
Capabilities of CODSTRAN Include determining the durability of composites with 
and without defects, determining structural responses due to mechanical and 
thermal loads, accurate prediction of stress states near defects (stress con- 
centrations), and prediction of ply and laminate level failure and fracture. 

The modules comprising CODSIRAN are: (1) the executive module, containing 

communication links to all other modules; (2) the I/O module; (3) the Analysis 
module; (4) the Composite Mechanics Module; and (5) the Fracture Mechanics 
module. 

The Analysis module Is the general purpose finite element computer code 
NAS1RAN and Is used to calculate both near- field and far-fleld stresses In a 
finite element model of the structural component or specimen. The Composite 
Mechanics (ref. 3) module generates laminate properties from constituent prop- 
erties (composite micromechanics) and uses Intraply failure and Interlaminar 
delamination criteria to check ply and Interply failure, respectively. The 
Fracture Mechanics module Is able to account for both ply level fracture and 
laminate fracture. The modified distortion energy principle and/or a general 
quadratic surface fit are used to Indicate combined ply level fracture. Lami- 
nate level failure Is assumed to occur when all the plies In the laminate have 
failed In all the failure modes. 

To computationally simulate progressive fracture, CODSTRAN uses an Intera- 
tlve procedure whereby a load Is applied Incrementally to the finite element 
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model of the structure being modeled. The response of the structure to the 
load can be no damage, damage, or destruction of an element(s) (local frac- 
ture). Based upon this response, the load Increment Is updated as follows: 

(1) If no damage Is predicted, the load Is updated by some predetermined load 
Increment; (2) If elements are damaged or destroyed (defect growth or local 
fracture), the same load Is reapplied with reduced material properties assigned 
to the damaged elements. Destroyed elements are purged from the finite element 
mesh, effectively defining progressive fracture. This load Is maintained until 
equilibrium In the structure Is achieved. Equilibrium, Is defined as the point 
where the structure, with Its updated geometry and modified material proper- 
ties, can sustain the applied load without the occurrence of further damage. 
This Iterative procedure Is continued through load Increments until global 
fracture of the structure occurs. Typical computationally simulated progres- 
sive fracture results are shown In figure 6. These results, perhaps more than 
any other, demonstrate the effectiveness of composite mechanical computational 
simulation. 


Fiber Composite Sandwich Thermostructural Behavior 

The use of fiber composites In space applications Is Increasing In a vari- 
ety of structural configurations. Sandwich structural configurations provide 
an effective application of fiber composites. The faces of the sandwich resist 
loads by membrane action which Is the most structural-efficient use of fiber 
composite thin laminates. In addition, composite sandwiches can be tailored 
to meet very close thermal distortion tolerances such as those required for 
communication satellite antennas. 

Designs to meet close thermal distortion tolerances are determined by 
detailed heat transfer, thermal, and structural analyses. These analyses 
require respective thermal and mechanical properties of the composite sandwich 
as well as temperature and moisture effects on these properties. Thermal prop- 
erties Include heat capacity, thermal conductivities and thermal expansion 
coefficients In the plane and through the thickness of the sandwich. The cor- 
responding mechanical properties Include normal and shear moduli and Poisson's 
ratios. Thermomechanical and corresponding thermostructural properties of 
sandwiches can In principle be measured. This Is practical only for one or at 
most a few sandwich configurations, limited temperature/moisture conditions and 
generally for selected properties. The alternative Is to computationally simu- 
late the thermomechanical behavior of composite sandwiches so that all the 
properties required for thermal and mechanical/structural analyses can be pre- 
dicted Including hygrothermal effects. Though this Is seemingly "straight- 
forward," It Is not available In the literature to the authors' knowledge. 
Recent research at NASA Lewis focused on developing computational methods for 
simulating the thermomechanical behavior of composite sandwiches (ref. 10). 
These methods use analyses with several levels of progressive sophistication/ 
simplification In conjunction with composite hygrothermomechanlcal theory. 

The several levels of sophistication Include: (1) three-dimensional 

finite element modeling of the honeycomb, the adhesive and the composite faces 
(fig- 7); (2) three-dimensional finite element modeling of the honeycomb assum- 
ing an equivalent homogeneous medium, the adhesive and the composite faces 
(fig. 8); (3) laminate theory simulation where the honeycomb (metal or compos- 
ite) Is assumed to be of plies with equivalent properties (fig. 8 without the 
finite element mesh); (4) derivation of approximate, simplified equations for 
simulating the honeycomb thermal and mechanical properties with an equivalent 
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homogeneous medium; and (5) hygrothermal effects from cryogenic to 250 °F with 
and without moisture. 

Typical results obtained from the several levels of sophistication/ 
simplification are summarized in table I. These results illustrate yet another 
application of computational simulation: an effective means to aid the devel- 

opment of approximate simple and explicit methods for predicting complex 
thermostructural behavior. 

High Temperature Metal Matrix Composite ( HT-MMC) Structure s 

The use of advanced high temperature composite materials in aerospace pro- 
pulsion systsems is becoming more attractive. The potential benefit is derived 
from the ability of these materials to allow higher operating temperatures 
resulting in Improved engine performance/efficiency and, at the same time, to 
enhance or at least maintain component durability. For example, qualitative 
investigations have identified the tungsten-f iber-reinforced superalloys (IFRS) 
as a potential first-generation composite aircraft turbine and Space Shuttle 
Main Engine (SSME) turbopump blade/vane material by virtue of the excellent 
combination of complementary properties provided by this family of advanced 
high-temperature composites. In order to adequately investigate the feasibil- 
ity of using TFRS for turbine engine and SSME components, computational simula- 
tion methods are needed which account for the complex physics associated with 
these applications including: nonlinear (stress-temperature-time dependent) 

and anisotropic constituent (f iber/matrix) material behavior, micro- and macro- 
heterogenity, and fiber/matrix Interdiffusion. 

Research activities at NASA Lewis over the past 5 yr lead to the develop 
ment of Nonlinear C0BS1RAN (Composite Blade Structural Analysis) for the com- 
putational simulation of HT-MMC blades (refs. 11 to 13). Nonlinear COBSTRAN 
embodies an upward Integration/downward decomposition strategy to predict 
thermomechanical structural response. This strategy comprises: (1) thermo- 

mechanical diffusion model, (2) nonlinear material constitutive models, (3) 
composite mechanics, and (4) global finite element structural analysis 
(fig. 9). Nonlinear COBSIRAN has been used to trace the behavior of HT-MMC 
components from fabrication through one or more missions. Results demonstrate 
that nonlinear COBSIRAN is computationally effective. It provides a direct 
means for predicting the complex behavior at all levels of the composite simu- 
lation; that is, from the constituent materials, through the several levels of 
composite mechanics, to the local/global structural response of complex HT-MMC 
structural components sub- ject to representative engine mission loads 
(fig. 10). 

Typical results obtained for the blade and mission shown in figure 10 are 
presented in figure 11. The Important observation is that it would be prac- 
tically impossible to determine this type of detailed behavior without computa- 
tional simulation. 


CONCLUSIONS 

Over the past 15 yr NASA Lewis Research Center has been developing spec- 
ialty methods for the analysis/design of composite structures. By virtue of 
their complexity, these methods are computational intensive and, by necessity, 
are in the form of Independent computer codes or are embedded as modules in 
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more general structural analysis codes. Collectively, the development of these 
various specialty methods constitute a new discipline Identified as computa- 
tional composite mechanics. These methods encompass all aspects of composite 
mechanics, Impact, progressive fracture, and component specific composite 
structural analysis. The application of computational simulation to select 
examples are briefly described In terms of (1) the Importance of the behavior 
simulated, (2) the simulation procedure, (3) typical results obtained, and (4) 
Interpretation and significance of these results. Collectively, these examples 
demonstrate the effectiveness, and perhaps the only coupled/integrated means of 
computational simulation for determining the complex structural response of 
composite structures In aerospace propulsion service environments. 
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TABLE I. - SUMMARY OF RESULTS FROM THREE DIFFERENT 
SIMULATION LEVELS OF PROGRESSIVE SIMPLIFICATION 


Property description 

Three-dimensional 
finite element 

Laminate 

theory 

Property 

Units 

Symbol 

Direct. 

Detailed 

model 

Homogen. 
core model 

Normal 

mpsl 

E 

XX 

1 

.0 

1 

.0 

1 .0 

modulus 

mpsl 


yy 

1 

.0 

1 

.0 

1 .0 


mpsl 


ZZ 


.09 


.09 

.09 

Shear 

mpsl 

G 

xy 


.39 


.38 

.38 

modulus 



yx 







mpsl 

G 

yz 


.071 


.071 

.071 




zy 







mpsl 

G 

zx 


.075 


.074 

.076 




xz 






Poisson 1 s 


V 

xy 


.396 


.315 

.314 

ratio 



yx 


.348 


.315 

.314 



V 

yz 


.121 


.043 





zy 


.073 


.005 




u 

zx 


.030 


.004 





xz 


.118 


.045 


Thermal 

y 1n./1n./°F 

a 

XX 

1 

.31 

1 

.21 

1 .21 

exp. coeff. 



yy 

1 

.38 

1 

.21 

1 .21 




ZZ 

— 

— 

14 

.6 

13.6 

Thermal heat 

Btu-1n./ 

K 

XX 


.191 


.191 

.190 


hr/in. 2/°F 


yy 


.181 


.183 

.182 




ZZ 


.043 


.043 

.041 
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3- ms FRAMES 


5-ms FRAMES 



KEVLAR/EPOXY 



T-300 EPOXY 



S-GLASS/EPOXY 

FIGURE 3. - DISPLACEMENT PROPAGATION DUE TO IMPACT. 


• DETERMINE REQUISITE PROPERTIES AT DESIRED CONDITIONS USING COMPOSITE 
MICROMECHANICS 

• RUN 3-D FINITE ELEMENT ANALYSIS ON ENF (MMF! FOR AN ARBITRARY LOAD 

• SCALE LOAD TO MATCH INTERLAMINAR SHEAR STRESS AT ELEMENT NEXT TO 
CRACK- TIP 

• WITH SCALED LOAD EXTEND CRACK AND PLOT STRAIN ENERGY RELEASE RATE 
VERSUS CRACK LENGTH 


• SELECT CRITICAL "G"l 
AND CRITICAL "a" J 

EXTENDED CRACK LENGTH, a 



• METHOD HAS VERSATILITY/GENERALITY 

FIGURE 4. - GENERAL PROCEDURE FOR SIMULATING INTERLAMINAR FRACTURE 
TOUGHNESS (END-NOTCH OF MIXED-MODE FRACTURE). 
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FIGURE 5. - SCHEMATICS OF TEST SPECIMENS FOR INTERLAMINAR FRACTURE 
TOUGHNESS. 



FIGURE 6. - COMPUTATIONALLY SIMULATED PROGRESSIVE FRACTURE 
AT VARIOUS PERCENTAGES OF THE FRACTURE LOAD (l5 402 lb; 

G r /E 1 0/+ 30/0/- 30/0 1 2 S ) . 
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pf POOk QUALITY 


FIGURE 7. - 3-D DETAILED FINITE ELEMENT MODEL OF COMPOSITE SAND- 
WHICH WITH HONEYCOMB CORE. 




FIGURE 8. - 3-D FINITE ELEMENT MODEL OF COMPOSITE SANDWICH WITH 
EQUIVALENT HOMOGENEOUS-MEDIUM CORE. 
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TO 

GLOBAL 

STRUCTURAL 

ANALYSIS 


/ 


COMPONENT 





FROM 

GLOBAL 

STRUCTURAL 

ANALYSIS 


LAMINATE 



LAMINATE 

THEORY 



LAMINATE I 
THEORY I 



COMPOSITE 

MICROMECHANICS 

THEORY 


\ 

CONSTITUENTS 




NONLINEAR > COMPOSITE 

MATERIAL r MICROMECHANICS 

MODEL | V THEORY 


o 


MATERIAL PROPERTIES 
P ■ P(o, T, t) 


FIGURE 9. - NONLINEAR COBSTRAIN SIMULATION OF HIGH TEMPERATURE METAL MATRIX 
COMPOSITES. 
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